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QCD at non-zero temperature: bulk properties and heavy quarks 
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Abstract: I review recent progress in lattice QCD at non-zero temperature wliith em- 
phasis on the calculations of equation of state and the properties of heavy quar anti-quark 
pairs at high temperatures. I also briefly discuss the deconfinement and chiral symmetry 
restoring aspects of the QCD transition at finite temperature 

O 1 Introduction 

^! It is expected that strongly interacting matter undergoes a transition in some temperature 

I interval from hadron gas to deconfined state also called the quark gluon plasma (QGP) 

Q> ' [1]. Creating deconfined medium in a laboratory is the subject of the large experimental 

program at RHIC [2] and is going to be the goal of the future heavy-ion program at LHC 
[3]. Attempts to study QCD thermodynamics on the lattice go back to the early 80's when 
c3 ! lattice calculations in SU (2) gauge theory provided the first rigorous theoretical evidence 

for deconfinement [U O [6]. The problem of calculating thermodynamic observables in 
D ' pure gluonic theory was solved in 1996 |7], while calculation involving dynamical quarks 

r^, ■ were limited to large quark masses and had no control over discretization errors [H [U [10] 

(see Refs. fH] W2\ for reviews). Lattice calculations of QCD thermodynamics with 
light dynamical quarks remained challenging until recently. During the past 5 years 
lO ! calculations with light u, d quarks have been performed using improved staggered fermion 

; actions [ia[Tan5l[l6l[I71tt8l[T9l|2niEIl[2a|23l[2l] (see Refs. [261 [ST].) 

\a ' To get reliable predictions from lattice QCD the lattice spacing a should be suffi- 

Q\ ■ ciently small relative to the typical QCD scale, i.e. Aqcdo. ^ 1. For staggered fermions, 

O I which are used for calculations at non-zero temperature, discretization errors go like 

C((aAQ(7£))^) but discretization errors due to flavor symmetry breaking turn out to be 
numerically quite large. To reduce these errors one has to use improved staggered fermion 
actions with so-called fat links [30]. At high temperature the dominant discretization er- 
rors go like {aTY and therefore could be very large. Thus it is mandatory to use improved 
discretization schemes, which improve the quark dispersion relation and eliminate these 
discretization errors. Lattice fermion actions used in numerical calculations typically im- 
plement some version of fat links as well as improvement of quark dispersion relation and 
are referred to as p4, asqtad, HISQ and stout. In lattice calculations the temperature 
is varied by varying the lattice spacing at fixed value of the temporal extent A^,-. The 
temperature T is related to lattice spacing and temporal extent, T = l/(A^a). There- 
fore taking the continuum limit corresponds to A^ — )■ cxd at the fixed physical volume. 
For the same reason discretization errors in the hadronic phase could be large when the 
temperature T is small. 

In this paper I am going to discuss lattice QCD calculation of the transition temper- 
ature, equation of state as well as different spatial and temporal correlation functions 
relevant for phenomenology of heavy quarks. 
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2 Equation of State 



The equation of state has been calculated with p4 and asqtad action on lattices with 
temporal extent iV^ = 4, 6 and 8 [TTl [HI [20] . In these calculations the strange quark mass 
was fixed to its physical value, while the light {u, d) quark masses 10 times smaller than 
the strange quark mass have been used. These correspond to pion masses of (220 — 260) 
MeV. The calculation of thermodynamic observables proceeds through the calculation 
of the trace of the energy momentum tensor e — 3p also known as trace anomaly or 
interaction measure. This is due to the fact that this quantity can be expressed in terms 
of expectation values of local gluonic and fermionic operators. The explicit expression 
for e — 3p in terms of these operators for p4 and asqtad actions can be found in Ref. 
[2U] . Different thermodynamic observables can be obtained from the interaction measure 
through integration. The pressure can be written as 



The lower integration limit To is chosen such that the pressure is exponentially small 
there. Furthermore, the entropy density can be written as s = (e + p)/T. Since the 
interaction measure is the basic thermodynamic observable in the lattice calculations it 
is worth discussing its properties more in detail. In Fig. [T]I show the interaction measure 
for p4 and asqtad actions for two different lattice spacings corresponding to = 6 and 
8. In the high temperature region, T > 250 MeV results obtained with two different 
lattice spacings and two different actions agree quite well with each other. Furthermore, 
recent calculations with HISQ actions also give results for the trace anomaly which are 
consistent with these [21]. Discretization errors are visible in the temperature region, 
where e — 3p is close to its maximum as well as in the low temperature region. At low 
temperatures the lattice data have been compared with the hadron resonance gas (HRG). 
As one can see the lattice data fall below the resonance gas value. This is partly due to 
the fact that the light quark masses are still about two times larger than the physical 
value as well as to discretization errors. The 0{{aAQCDy) discretization errors in the 
hadron spectrum are suppressed at high temperatures as the lattice spacing a is small 
there. Also hadrons are not the relevant degrees of freedom in this temperature region. 
But at small temperatures, where hadrons are the relevant degrees of freedom, these 
discretization effects are significant. It turns out, however, that the HRG model that 
takes into account the quark mass dependence and discretization errors in the hadron 
spectrum can describe the lattice data quite well [311 [32]. The large discretization errors 
in the low temperature region is the reason for the discrepancies with the stout results at 
temperatures T < 200 MeV [23]. The lattice results for the trace anomaly obtained with 
stout action are also different in the high temperature limit [22]. Namely, (e — 3p)/T^ 
calculatied with stout action is about 50% smaller in the high temperature region. To 
resolve this problem calculation at finer lattice spacing are needed. 

The pressure, the energy density and the entropy density are shown in Fig. [2l The 
energy density shows a rapid rise in the temperature region (185 — 195) MeV and quickly 
approaches about 90% of the ideal gas value. The pressure rises less rapidly but at the 
highest temperature it is also only about 15% below the ideal gas value. In the previous 
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Figure 1: The interaction measure calculated with pA and asqtad actions in entire 
temperature range (left) and at low temperatures (right) from Ref. [20]. The dashed and 
dashed-dotted lines are the prediction of hadron resonance gas (HRG) with all resonances 
included up 2.5GeV (dashed) and l.SGeV (dashed-dotted), respectively. 

calculations with the p4 action it was found that the pressure and energy density are 
below the ideal gas value by about 25% at high temperatures [TU]. Possible reason for 
this larger deviation could be the fact that the quark masses used in this calculation 
were fixed in units of temperature instead being tuned to give constant meson masses 
as lattice spacing is decreased. As discussed in Ref. [JT] this could reduce the pressure 
by 10 — 15% at high temperatures. In Fig. [2] I also show the entropy density divided 
by the corresponding ideal gas value and compare the results of lattice calculations with 
resummed perturbative calculation [331 |3l] well as with the predictions from AdS/CFT 
correspondence for the strongly coupled regime [35] . The later is considerably below the 
lattice results. Note that pressure, energy density and the trace anomaly have also been 
recently discussed in the framework of resummed perturbative calculations which seem 
to agree with lattice data quite well [36j . 

The differences between the stout action and the p4 and asqtad actions for the trace 
anomaly translates into the differences in the pressure and the energy density. In par- 
ticular, the energy density is about 20% below the ideal gas limit for the stout action. 



3 Fluctuations of conserved charges 

The pressure at non-zero chemical potentials can be evaluated using Taylor expansion. 
The Taylor expansion can be set up in terms of quark chemical potentials or in terms of 
chemical potentials corresponding to baryon number B, electric charge Q and strangeness 
5* of hadrons. The expansion coefficients in quark chemical potential xUds 
hadronic ones Xsqs s-^e related to each other. For a more detailed discussion on this 
see Ref. [37]. Taylor expansion can be used to study the physics at non-zero baryon 
density. However, the expansion coefficients also give important information about fluc- 
tuation of conserved charges, which turn out to be sensitive probes of deconflnement 
and chiral symmetry restoration. Fluctuations of conserved charges have been studied 
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Figure 2: The energy density and the pressure as function of the temperature (left), and 
the entropy density divided by the corresponding ideal gas value (right). The dashed 
lines in the right panel correspond to the resummed perturbative calculations while the 
solid black line is the AdS/CFT result. 



in detail using asqtad and p4 actions [I3l [20l |37]. In Fig. |3] I show the quadratic and 
quartic fluctuations of B, Q and S. Fluctuations are suppressed at low temperatures 
because conserved charges are carried by massive hadrons, but rapidly grow in the tran- 
sition region T = (185 — 195) MeV as consequence of deconfinement. At temperatures 
T > 300 MeV fluctuations are reasonably well described by ideal gas of quarks. Quartic 
fluctuations exhibit a peak near the transition region, which at sufficiently small quark 
masses can be related to the chiral transition and the corresponding critical exponents 
(see discussion in Refs. [20]). As one can see from the figure quadratic fluctuations show 
visible cutoff dependence at low temperatures, while the cutoff dependence is small for 
T > 250 MeV. 

The transition from hadronic to quark degrees of freedom can be particularly well 
seen in the temperature dependence of the Kurtosis, which is the ratio of quartic to 
quadratic fluctuations. This quantity can be also measured experimentally. In Fig. H] 
I show the Kurtosis of the baryon number as a function of the temperature. At low 
temperatures it is temperature independent and is close to unity in agreement with the 
prediction of the hadron resonance gas model. In the transition region it sharply drops 
from the hadron resonance gas value to the value corresponding to an ideal gas of quarks. 
Since the fluctuations of conserved charges are so well described by ideal quark gas it is 
interesting to compare the lattice results for quark number fluctuations with resummed 
perturbation theory. Quadratic quark number fluctuations are also called quark number 
susceptibilities and have been calculated in resummed perturbation theory |3Sl ESI SD] • 
The comparison of the resummed perturbative results with lattice calculations performed 
with p4 action [12] for the quark number susceptibility divided by the ideal gas value is 
shown in Fig. HJ The figure shows a very good agreement between the lattice calculations 
and resummed perturbative results. It is also interesting to compare findings of the 
lattice calculations with the results obtained in strongly coupled gauge theories using 
AdS/CFT correspondence. The result from AdS/CFT calculations shown in the figure 
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as the solid black line is significantly below the lattice results 0. Let me finally note that 
the discretization errors in the lattice calculations of the fluctuations are small at high 
temperatures when calculated with the pA action. The same is true for asqtad action 
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Figure 3: The quadratic (left) and quartic (right) fluctuations of B, Q and S as function 
of the temperature [37]. The open symbols correspond to the results obtained on A^.^ = 4 
lattices, while the flUed symbols to the ones obtained on A^,- = 6 lattices. 
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Figure 4: The Kurtosis of baryon number (left) and the quark susceptibility compared 
with resummed perturbatlve calculations (right). The open symbols on the right plot 
correspond to asqtad results for strange quark number susceptibility [I3] . The resummed 
perturbatlve results were obtained in next-to-leading approximation (NLA) [39]. 



4 Chiral and deconfinement transition 

The flnite temperature transition in QCD has aspects related to deconflnement and chi- 
ral symmetry restoration. Deconflnement aspects of the transition are related to color 
screening and, for inflnitely heavy quarks, also to the center symmetry. The order param- 
eter for deconflnement is the Polyakov loop which is related to the free energy of isolated 
static quark. The breaking of the center symmetry is signaled by non-zero value of the 
Polyakov loop. In the opposite limit of massless quarks QCD has the chiral symmetry. 

^the conserved charges considered in these calculations are not exactly the quark numbers, see dis- 
cussion in Ref. [26] . 
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Figure 5: The subtracted chiral condensate (left) and the renormahzed Polyakov loop 
(right) calculated for different actions [121 Ell [21] • 

The quark condensate (ipip) is the order parameter for this symmetry. The chiral sym- 
metry is broken in the vacuum and expected to be restored at high temperatures. In Fig. 
|5]the renormahzed Polyakov loop and the subtracted chiral condensate are shown. Both 
quantities show a smooth change in the transition region which is consistent with the 
fact that the finite temperature transition is an analytic crossover and not a true phase 
transition [13]. It is not clear to what extent the center symmetry is relevant for the 
physical values of the quark masses and it is difficult to define a deconfinement transition 
temperature. The universal aspects of the chiral transition seem to be relevant for the 
range of quark masses which include the physical light quark mass and the corresponding 
transition temperature can be defined. Calculations with the p4 action on A^,- = 4 and 
6 lattices gave an estimate Tc = 192(4) (7)MeV for the chiral transition temperature |15] 
which is significantly larger than the value of about 155MeV obtained for stout action. 
This is due to large flavor symmetry breaking for the p4 action. As one can see from Fig. 
|5]the discrepancies between different actions are reduced when considering larger Nr or 
the HISQ action where the effects of flavor symmetry breaking are much smaller. It is 
expected that all actions will give the same result in the continuum limit resulting in a 
transition temperature of (150 — 160)MeV. Discrertization effects are also visible for the 
renormahzed Polyakov loop in the low temperature region for p4 and asqtad calculations 
with Nt- = 8. However, the A^^ = 12 asqtad, HISQ and stout data are all agree. Thus, 
discretization errors are under control for this quantity. 

5 Correlation functions of static quark anti-quark 
pair 

One of the most prominent feature of the quark gluon plasma is the presence of chro- 
moelectric (Debye) screening. The easiest way to study chromoelectric screening is to 
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Figure 6: The singlet free energy Fi{r, T) calculated in Coulomb gauge on 16^ x 4 lattices 
(left) and the combination Fi{r,T) — Foo{T) as function of rT (right). The solid black 
line is the parametrization of the zero temperature potential. 

calculate the singlet free energy of static quark anti-quark pair (for recent reviews on this 
see Ref. [Hj SS]), which is expressed in terms of correlation function of temporal Wilson 
lines in Coulomb gauge 

exp(-Fi (r,T)/T) = ^Ti {W{r)W^m. (2) 

L = TrW is the Polyakov loop. Instead of using the Coulomb gauge the singlet free 
energy can be defined in gauge invariant manner by inserting a spatial gauge connection 
between the two Wilson lines. Using such definition the singlet free energy has been 
calculated in SU{2) gauge theory 06]. It has been found that the singlet free energy 
calculated this way is close to the result obtained in Coulomb gauge [16]. The singlet 
free energy turned out to be useful to study quarkonia binding at high temperatures 
in potential models (see e.g. Refs. [HI IHl [501 ED E2])- The singlet free energy also 
appears naturally in the perturbative calculations of the Polyakov loop correlators at 
short distances [47]. 

The singlet free energy was recently calculated in QCD with one strange quark and 
two light quarks with masses corresponding to pion mass of 220MeV on 16^ x 4 lattices 
|53] . The numerical results are shown in Fig. [61 At short distances the singlet free 
energy is temperature independent and coincides with the zero temperature potential. 
In purely gluonic theory the free energy grows linearly with the separation between the 
heavy quark and anti-quark in the confined phase. In presence of dynamical quarks the 
free energy is saturated at some finite value at distances of about 1 fm due to string 
breaking |1H [SH [57]. This is also seen in Fig. [6] Above the deconfinement temperature 
the singlet free energy is exponentially screened at sufficiently large distances \5E\ 
with the screening mass proportional to the temperature , i.e. 

Fi(r,T) = Foo(T) - l^!i^exp(-m,5(T)r), ~ T. (3) 
6 Airr 

Therefore in Fig. [6] we also show the combination Fi{r,T) — Foo(T) as a function of rT. 
As one can see from the figure this function shows an exponential fall-off at distances 
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rT > 0.8. The fact that the slope is the same for all temperatures means that ttid ~ T, 
as expected. 



6 Temporal meson correlars and spectral functions 

Information on hadron properties at finite temperature as well as the transport coef- 
ficients are encoded in different spectral functions. In particular the fate of different 
quarkonium states in QGP can studied by calculating the corresponding quarkonium 
spectral functions (for a recent review see Ref. [15]). On the lattice we can calculate cor- 
relation function in Euclidean time. This is related to the spectral function via integral 
relation 

nf ^\ ( ^\ i^i ^\ cosh(u;(r- 1/2T)) 

Given the data on the Euclidean meson correlator Gir^ T) the meson spectral function 
can be calculated using the Maximum Entropy Method (MEM) [58]. For charmonium 
this was done by using correlators calculated on isotropic lattices |59l [60] as well as 
anisotropic lattices [611 [621 [63] in the quenched approximation. It has been found that 
quarkonium correlation function in Euclidean time show only very small temperature 
dependence [60l [63] . In other channels, namely the vector, scalar and axial- vector chan- 
nels stronger temperature dependence was found [Sni[n3|- The spectral functions in the 
pseudo-scalar and vector channels reconstructed from MEM show peak structures which 
may be interpreted as a ground state peak [601 ED [62]. Together with the weak temper- 
ature dependence of the correlation functions this was taken as strong indication that 
the IS charmonia {r]c and J/ip) survive in the deconfined phase to temperatures as high 
as 1.6Tc [601 [nH [62] . A detailed study of the systematic effects show, however, that the 
reconstruction of the charmonium spectral function is not reliable at high temperatures 
[63] . in particular the presence of peaks corresponding to bound states cannot be reli- 
ably established. Presence of large cutoff effects at high frequencies also complicates the 
analysis [66]. The only statement that can be made is that the spectral function does 
not show significant changes within the errors of the calculations. Recently quarkonium 
spectral functions have been studied using potential models and lattice data for the sin- 
glet free energy of static quark anti-quark pair [501 ED ES]. These calculations show 
that all charmonium states are dissolved at temperatures smaller than 1.2Tc, but the 
Euclidean correlators do not show significant changes and are in fairly good agreement 
with available lattice data both for charmonium [601 ES] and bottomonium [63l [M] . This 
is due to the fact that even in absence of bound states quarkonium spectral functions 
show significant enhancement in the threshold region [49] . Therefore previous statements 
about quarkonia survival at high temperatures have to be revisited. Exploratory calcu- 
lations of the charmonium correlators and spectral functions in 2-flavor QCD have been 
reported in Ref. [65] and the qualitative behavior of the correlation functions was found 
to be similar. 

The large enhancement of the quarkonium correlators above deconfinement in the 
scalar and axial-vector channel can be understood in terms of the zero mode contribution 
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[i9| [67] and not due to the dissolution of the IP states as previously thought. Similar, 
though smaller in magnitude, enhancement of quarkonium correlators due to zero mode 
is seen also in the vector channel Here it is related to heavy quark transport 

[18| [68] . Due to the heavy quark mass the Euclidean correlators for heavy quarkonium 
can be decomposed into a high and low energy part G{t,T) = Giow{t,T) + Ghigh(T, T) 
The area under the peak in the spectral functions at zero energy u giving the 
zero mode contribution to the Euclidean correlator is proportional to some susceptibility, 
GI^^{t,T) ~ Tx^T), which have been calculated on the lattice in Ref. pH]- It is 
natural to ask whether the generalized susceptibilities can be described by a quasi-particle 
model. The generalized susceptibilities have been calculated in Ref. [70] in the free 
theory. Replacing the bare quark mass entering in the expression of the generalized 
susceptibilities by an effective temperature dependent masses one can describe the zero 
mode contribution very well in all channels [SS]- 

While temporal correlators are not sensitive to the change in the spectral functions 
spatial quarkonium correlation functions could be more sensitive to this. Recent lattice 
calculations show indication for significant change in spatial charmonium correlators 
above deconfinement [TT] . 

The spectral function for light mesons has been calculated on the lattice in quenched 
approximation [72 | [71 ] [75]. However, unlike in the quarkonia case the systematic errors 
in these calculations are not well understood. 

7 Conclusions 

In recent years significant progress has been achieved in studying strongly interacting 
matter at high temperatures using lattice QCD. Equation of state and transition tem- 
perature have been calculated at several lattice spacings allowing for controlled continuum 
extrapolations. Current lattice data suggest a transition temperature for chiral symmetry 
restoration of (150 — 160) MeV. Fluctuations of conserved charges and spatial correlation 
functions provide detailed picture of the high temperature medium, namely the nature 
of the relevant degrees of freedom (partonic or hadronic) and color screening. Temporal 
meson correlation functions have been also studied in lattice QCD but the numerical data 
are not precise enough to provide detailed information of the corresponding meson spec- 
tral functions. However, these lattice data are useful to constrain the model calculations 
of the meson spectral functions. 
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